Single-Species Three-Particle Reactions in One Dimension 
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Renormalization group calculations for fluctuation-dominated reaction-diffusion systems are gen- 
erally in agreement with simulations and exact solutions. However, simulations of the single-species 
reactions 3A — > (0, A, 2A) at their upper critical dimension d c = 1 have found asymptotic densities 
argued to be inconsistent with renormalization group predictions. We show that this discrepancy 
is resolved by inclusion of the leading corrections to scaling, which we derive explicitly and show to 
be universal, a property not shared by the A + A — > (0, A) reactions. Finally, we demonstrate that 
two previous Smoluchowski approaches to this problem reduce, with various corrections, to a single 
theory which surprisingly yields the same asymptotic density as the renormalization group. 



I. INTRODUCTION 



Reaction-diffusion systems are known to be strongly 
dependent on fluctuations when the spatial dimension 
d is at or below an upper critical dimension d c . This 
fluctuation-dominated case has been treated by field- 
theoretic renormalization group (RG) methods for a wide 
variety of reaction types and conditions, as recently re- 
viewed Q. Comparison of the RG results with exact 
solutions and simulations has generally yielded agree- 
ment or at least consistency, as detailed in examples 
given below. However, simulations of the single-species 
3A — > (0, A, 2A) reactions at the upper critical dimen- 
sion a = a r appear to be inconsistent with 
RG predictions Q . This discrepancy is noteworthy since 
these reactions present one of the simplest and directly 
testable cases. In the current work we demonstrate that 
there is no discrepancy. 

The general single-species reaction-diffusion decay 
kA — > IA with integers I < k has an upper critical dimen- 
sion d c = 2/(k — 1). Above the upper critical dimension 
the density n(t) follows the rate equation, dtn = —\n k , 
which gives the asymptotic decay n ~ f -1 /^ -1 ) with 
an amplitude that depends on the nonuniversal rate con- 
stant A. Below the upper critical dimension particle anti- 
correlations neglected by the rate equation become rel- 
evant, giving rise to a slower decay n ~ Ak,e(Dt)~ d / 2 , 
with Ak,i a universal constant dependent only on the re- 
action type (k and £) and the dimension d. RG methods 
have been used to show the exponent to be exact [j| , and 
to demonstrate the universality of A^e and provide an 
expansion in powers of d c — d for the amplitude 0. 

The amplitude expansion may be tested by comparison 
with exact solutions and simulations. For example, solv- 
able realizations of the A + A — > model in one dimen- 
sion give the decay amplitude of A 2 ,o = l/v87r — 0.199 
0, . The RG expansion in e = d c — d = 1 gives instead 
A 2fi ^ 0.080 + 0.029 + . . . 0. Evidently the truncated 
perturbative RG is of little accuracy when e = 1. How- 
ever, if the diffusive transport is replaced by long-range 
hops, the upper critical dimension can be continuously 
lowered from two to one, allowing for a truly small e 
in spatial dimension d — 1 8]. In this case, the theo- 



retical amplitude compares well with simulations. Also 
noteworthy is that the ratio A 2 ,\IA 2 fl = 2 found from 
exact solutions |6j, |9j] is also an exact result from the RG 
calculation, i.e., a field rescaling transformation in the 
field theory shows this ratio to hold to all orders in the e 
expansion Thus simulations and exact solutions for 
d < d c are in general agreement with the field-theoretic 
RG approach. 

For the borderline case of d = d c , the density is pre- 
dicted to decay with the rate equation exponent, but with 
logarithmic corrections and a universal amplitude 
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The amplitude in this case is given explicitly, rather than 
perturbatively, as 

Ak ' e= { Mk—J) ) ■ (2) 

This explicit result provides then a strong test for the RG 
calculation. An exact solution is available for a particular 
realization of the A + A — > (0, A) reactions in d = d c = 2 
0, with values ^2,0 = 1/8tt and A 2 ,i = 1/Att that match 
the RG results, eq. 101. 

However, the predictions for the 3^4 — > IA reactions at 
the upper critical dimension d c = 1 have been the source 
of some controversy. Simulations have demonstrated log- 
arithmic corrections in d = 1, but with amplitudes that 
differ from the renormalization group predictions. Specif- 
ically, while eq. © gives A 3J ~ 0.21, 0.26, and 0.37 for 
I = 0, 1, and 2 respectively, simulations have reported 
values A 3:0 ~ 0.26 3j, A 3A ~ 0.76, and A^ 2 ^ 0.93 0. 
Further, the 3A — * simulations were found to be con- 
sistent with a version of Smoluchowski theory adapted 
for three-particle reactions That is, an approximate 
theory appears to agree better with the simulations than 
the RG calculation, which in principle involves no ap- 
proximations. 

To address this discrepancy, we revisit all of the field- 
theoretic approach, simulations, and Smoluchowski the- 
ory. Our main results are as follows. First, we demon- 
strate with RG methods that the leading corrections to 
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the asymptotic density are universal, a surprising prop- 
erty not shared by the A + A — ► (0, A) reactions at their 
upper critical dimension d c = 2. Our result is 
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with given by eq. J5J and B3 j computed in sectionlTTl 
below. Explicitly, 
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The next term in the expansion is nonuniversal. 

This universal leading correction is quite significant in 
the time range available to simulations, of the order of 
half the magnitude of the asymptotic density. Conse- 
quently, including this term makes the RG predictions 
consistent with the previous 3A — ► simulations but 
not with the 3 A — > (A, 2 A) simulations Q. 

This motivated us to conduct our own simulations, 
which are presented in section |TTT| Our simulation data 
is consistent with that of Ref. [3j for the 3A — > reac- 
tion, but not compatible with the data of Ref. Q for the 
3A — > (A, 2A) reactions, which we believe to be in err. 
Our simulation data is not capable of directly confirming 
the predicted amplitudes due to remaining, slow tran- 
sients. However, the data shows no discrepancy with the 
RG predictions for all three cases. 

In order to have a test of the RG predictions that does 
not contain slow transients, we turn in section HV1 to the 
field rescaling transformation in the field theory, which 
we use to predict relations between pure and mixed re- 
actions. Our simulations verify that these relations hold 
with high accuracy. 

Finally, we turn to the Smoluchowski approach, which 
is presented in section [V] We begin by demonstrating 
that the Smoluchowski solution for the — > (0, A) re- 
actions in d — 2 not only exhibits the logarithmic correc- 
tions, but also gives the correct density amplitude A2/. 
We believe this to be a new result. This is also germane 
to the 3A — > IA reaction in d = d c = 1 because both 
Smoluchowski approaches to this problem in the litera- 
ture 0, ITo| are constructed via a quasi-two-dimensional 
approach. We show that these two approaches, when 
various omitted factors are included, reduce to the same 
Smoluchowski theory, and further that this theory repro- 
duces the same asymptotic density as the RG approach. 

A summary is presented in section IVT1 



II. 



RENORMALIZATION GROUP 
CALCULATION 



Our presentation in this section will follow closely the 
formalism developed in The Doi-Peliti mapping of 
reaction-diffusion systems to a field theory is by now a 



standard technique Q, 0, ^3 1 gi ym g f° r the 3 A 
reaction the action 



S = / d d xdt 
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The first term in the integrand corresponds to the dif- 
fusion process and provides the propagator for the field 
theory. The higher order terms correspond to the reac- 
tion and provide vertices in the diagrammatic expansion. 
The Poisson initial conditions are reflected in the initial 
term no4>{t = 0). 

Here Ao is a nonuniversal, bare coupling constant asso- 
ciated with the microscopic reaction rate, while the coef- 
ficients ci and C2 depend only on £. These coefficients are 
determined by a field shift: the Doi-Pcliti mapping first 
gives the action in terms of fields 4> an d <f>, with the cou- 
pling terms Ao(</> 3 — (j> e )cf> 3 . To eliminate a "final term" 
4>(tf) that complicates the calculations (see 0,0 for de- 
tails) the field shift <j) = 1 + cf> is then employed, resulting 
in the coupling terms in eq. (JSJl, with the coefficients 



C2 = \{3-t)(2 + l). 



(0) 



The renormalization of this action Q, Q , necessary to 
obtain finite calculations for d < I, is relatively straight- 
forward, requiring only renormalization of the coupling 
constant. An arbitrary normalization time io is used to 
define a dimensionless, renormalized coupling constant 
qr. The renormalization group flow, for t 3> to, is de- 
scribed by the Callan-Symanzik equation 

n(t,n ,g R ,t ) = (to/t) d/2 n(t ,no(t),g R (t),t ) (7) 



with the running initial density 

no(t) - n a {t/t ) d/2 
and running coupling for d = d c = 1 given by 

9R(t) 



V3ln(t/r) 



(8) 



(9) 



where r is some nonuniversal time constant related, via 
the diffusion constant, to the short distance cutoff of the 
field theory The coefficient in eq. 10 is determined by 
the loop integrals in the definition of the renormalized 
coupling constant. The logarithmic time dependence is 
characteristic of marginal operators at the upper critical 
dimension. 

A loop expansion in terms of the bare coupling A and 
initial density no is used for the right hand side of eq. 
to give the asymptotic t ^ to density. For example, 
the tree-level (zero loop) diagrams may be summed by a 
Dyson equation 4], giving 
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The coupling Xq/D is converted to ga and then replaced 
by the running coupling c/r. The initial density n is 
replaced by ho, which grows as 

t d/2 

via eq. JBJ • This lat- 
ter step allows us to keep only the large no limit in the 
unrenormalized density calculation. Any corrections due 
to finite no will renormalize to subasymptotic contribu- 
tions to the density. Finally, the time is set to to and the 
prefactor in eq. J7J is included, giving the renormalized 
contribution 
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The higher order diagrams give renormalized contribu- 
tions at order (Dt) -1 / 2 ^^ , where m is the number 
of loops. Since ga ~ 1/lni for large t, these represent 
subleading terms to the asymptotic density. Thus with 
eq. © for c\ we obtain the leading order asymptotic den- 
sity reported in jj], and given by eq. (|T|). 

Up to this point we have summarized results from pj. 
Now we demonstrate that the leading corrections to the 
asymptotic density are themselves universal. The renor- 
malized loop expansion with gn ~ l/ln(£/r) takes the 
form 



(Dt) 1/2 n{t) = As/WtJT) + B + C/sfirttfF) + ... (12) 



where A, B, C, . . . are universal coefficients. Nonuniver- 
sal terms, apart from r above, are suppressed by negative 
powers of time. Reor gan izing the e xpa nsion in terms of 
lni via we observe 

that the nonuniversal r dependence does not enter until 
order 1 / Vhi t. Hence, B^e represents a universal leading 
correction to the asymptotic density. 

This is in contrast to the A + A — ► (0, A) reaction at 
the upper critical dimension, for which the renormalized 
loop expansion has the form (Dt)^ 1 [A ln(i/r) + B + . . .], 
and the r dependence enters the leading correction. 

Next we determine the coefficient. As presented in 
|4( , the summation of all one-loop diagrams requires the 
use of the tree-level density n^(t), given by eq. |JTDJ, 
and the tree- level response function (x, t\, t 2 ) = 
(<f>(x,t2)4>{0,ti))o, where the subscript on the angle 
brackets indicates a tree-level average. With the Fourier 
transform f(p) = J dxe~ lpx f(x) we obtain 



3/2 
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in the large no limit, which is sufficient for our purposes. 
From the one-loop diagram contribution is given by 



t rt 2 

(t) = 6 / dt 2 I dti 



dp 
2^ 



G^(0,t,i 2 )(- Cl Ao)n (0) (i2)G( )(p,t 2 ,t 1 )G( )(-p,t 2 ,t 1 )(~ C2 A )[n( )(i 1 )] 3 (14) 
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The factor of six results from the combinatorics of at- 
taching the response functions to the vertices. Evaluat- 
ing eq. i|14|) and using the result in the right-hand side of 
the CS equation Q gives the one-loop contribution 



n£>(t) 
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Substituting for ci and c 2 via eq. JBJl gives the value for 
B 3 £ reported in eq. Q- 



III. SIMULATIONS 

In order to test our density calculations, we simulate 
the 3A — > (0, A, 2A) reactions in one dimension. We 
employ synchronous dynamics, in which the particles are 
restricted to be on all even or all odd numbered sites at 
a given time. All particles are updated simultaneously, 
each particle randomly hoping one site left or right, which 
corresponds with a lattice spacing of unity to a diffusion 
constant D = 1/2. Subsequent reactions occur between 
particles on the same site until there are no more than 
two particles per site. This completes a single time step. 



For initial conditions, each lattice site is randomly as- 
signed 0, 1, or 2 particles with probabilities po, Pi, and 
P2 respectively. We choose the pi as follows: first, the 
particle numbers are determined by a Poisson distribu- 
tion with average no- Then these particles are allowed to 
react until there are no more than two particles per site, 
which gives the site occupation probabilities (po, Pi, P2) 
as functions of no and £. In order to reduce initial density 
transients, we then take the limit no — > 00, which gives 



(P0,Pl,P2) 




(16) 



These probabilities are used as the initial conditions for 
our simulations. All simulations were conducted on a 
lattice of size 2 23 ~ 8.4 x 10 6 for 10 8 time steps, with 10 
to 20 independent runs for each case. 

In Fig. H w e show the densities for the three reactions, 
plotted as n{t) 2 Dt versus In t. The data are clearly incon- 
sistent with the rate equation result, n ~ t" 1 / 2 , suggest- 
ing logarithmic corrections. According to eq. Q the den- 
sity plotted with these axes should approach a straight 
line with slope A 2 t While the data show seemingly little 
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FIG. 1: Plot of n(t) 2 Dt versus lnt for the 3A -> (0, A, 2A) 
reactions, as labeled. The three curves for the 3A — » re- 
action correspond, from top to bottom, to initial densities 
no — > oo, no =0.1 and no = 0.01. The rate equation predic- 
tion n ~ i _1//2 would correspond to a horizontal line. 
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FIG. 2: Plot of the density scaled by its expected time de- 
pendence versus time for the 3A — > reaction. The data is 
averaged over 20 independent runs. Also shown are the RG 
predictions for the asymptote, a constant on this plot, and 
the asymptote with universal leading corrections. 



curvature at late times, we present evidence below that 
suggests the asymptotic slopes have not been reached. 

Also illustrated in Fig. is the universality with re- 
spect to the initial density, shown for the 3A — > reac- 
tion. Initial densities of no = 0.1 and no = 0.01 produce 
only transient deviations from the no — * oo data. The 
duration of the transient grows with decreasing no- 

Next we show, in turn, the same data for each of the 
pure reactions, as compared to the RG predictions. In 
Fig. |21 the density is multiplied by y/Dt/hit, the in- 
verse of the expected time dependence, and this is plot- 
ted versus time. According to the RG prediction, the 
data should approach the constant value asymptot- 
ically. This value included on the plot, and is seen to be 
well below the data. Also plotted is the universal sum 
of the asymptotic density and leading corrections. It is 
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FIG. 3: Plot of the 3 A — > A simulation data, with axes and 
RG predictions as described in Fig. |5] The data is averaged 
over 15 independent runs. 
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FIG. 4: Plot of the 3^4 — > 2A simulation data, with axes and 
RG predictions as described in Fig. The data is averaged 
over 10 independent runs. 



clear the leading corrections play a significant role in the 
time range accessible to simulation. It is important to 
note that the upper RG curve has the lower constant as 
its asymptote. 

Our simulations appear to be consistent with those of 
Ref. Q • We further conclude that the simulations are not 
in conflict with the RG calculations. The slow, logarith- 
mic decay of the nonuniversal transient is responsible for 
the remaining discrepancy between the simulation data 
and the RG. Extending the simulations to asymptotia 
would appear to be impossible. 

Similar plots for the I = 1 and I — 2 reactions are 
shown in Figs. and 0] In each case, the data fall well 
above the RG asymptote, but below the RG asymptote 
with the universal leading corrections. Here we note, 
observing the vertical scale, that our data are not con- 
sistent with the amplitudes 0.76 (£ = 1) and 0.93 (£ — 2) 
reported in |2|. 
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IV. MIXED REACTIONS AND THE FIELD 
RESCALING TRANSFORMATION 

The RG field rescaling transformation, described be- 
low, predicts relationships between various mixed and 
pure reactions, similar to the predicted amplitude ratio 
for the k = 2 case. These predictions can be tested by 
simulations, and furthermore the tests are not plagued by 
the logarithmically slow transients found in the previous 
section. 

As discussed in Q , the action eq. 10 includes the case 
of mixed reactions, where different reactions may occur 
according to specified probabilities. Specifically, let qt be 
the probability that I particles remain after three parti- 
cles meet, with (70 + 91+92 = 1- These reaction rates 
carry linearly through the Doi-Pcliti formalism from the 
master equation to the field theory, giving the same ac- 
tion JSJ), with the coefficients 

2 

ci - ^(3-^ = l + 2g + ?i (17) 

2 

c-2 = ^i(3-£)(2 + £) g£ = 2 + go + <7i (18) 

e=o 

when expressed in terms of go & n( l Qi- We can write the 
universal density amplitudes in eq. J2J in terms of c\ and 
C2 as 

A ( ^ r _9V2^ 

M *-\*^) ' ^-"e^T' (19) 

via eqs. H10|) and l|15(l . Thus we have a prediction for 
the density for any mixed reaction as well as the pure 
reactions. 

A field rescaling transformation on the action eq. © 
may be used to predict relations between various pure 
and mixed reactions. The transformation <fi — > 6</> and 
<fi — > 6~ 1 </> leaves the diffusion part of the action un- 
changed, while transforming the coupling prefactors and 
initial density according to 

ci — > b~ 2 c\ C2 — > b~ 1 C2 no — > brio. (20) 
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FIG. 5: The density n(t) for the 3A — + 2A reaction, as well 
as for two mixed reactions as labeled (see text). Also plotted 
is the 3 A — > 2 A data scaled by factors of 3/4 and 2/3, which 
overlays the mixed reaction data as predicted by the field 
rescaling transformation. 

ultimately predicts that the density from different sys- 
tems may be closely related by the field rescaling trans- 
formation, even at times accessible to simulation. This 
provides an additional test of the RG formalism. 

In Fig. we show simulation data for the pure 
3A — ► 2A reaction along with two mixed reactions de- 
fined by the probabilities (<7o,9i,92) = (1/4,3/4,0) and 
(1/9,5/9, 1/3). According to the field rescaling transfor- 
mation, these mixed reactions are both related to the 
1 = 2 case with rescaling factors 6 = 2/3 and 6 = 3/4 
respectively. The pure reaction data is also plotted mul- 
tiplied by 6, which is seen to match quite well with both 
cases of mixed reaction data. Evidently the nonuniversal 
time constant r, which would be the only source of a slow 
transient difference, is essentially the same in the various 
reactions. We note that the reactions were implemented 
with the same microscopic rules, i.e., synchronous dy- 
namics with a maximum occupancy of two particles per 
site and reaction occurring only between particles on the 
same site. 



V. SMOLUCHOWSKI THEORY 



Thus, for example, the pure £ = 2 action, with c\ = 1 
and C2 = 2, can be transformed by a factor 6 = 2/3 
to a mixed reaction with (qo,<7i,<72) = (1/4, 3/4, 0). As 
a result, the field-theory loop expansion for density in 
the 3A — > 2A reaction, when multiplied by 6, will match 
that of the mixed reaction to all orders in the no — > oo 
limit. The only difference in the resulting renormalized 
densities will be a short-lived transients from the initial 
densities and irrelevant couplings, and a longer-lived but 
relatively small transient from differing nonuniversal time 
constants r in the marginal renormalized coupling. Since 
the r values may be quite similar or identical for mod- 
els which implement the reaction the same way, the RG 



Smoluchowski theory is a mean-field theory based on 
the correlation function, or conditional density, as com- 
pared to rate equations which are based on the density. 
Interestingly, Smoluchowski theory is capable of captur- 
ing much of the behavior of the d < d c , fluctuation- 
dominated systems missed by the rate equations: for the 
A + A — > (0, A) reactions it gives the correct upper crit- 
ical dimension d c = 2, correct decay exponents above 
and below d c , and even exhibits logarithmic corrections 
at the critical dimension. However, the approach does 
involve an uncontrolled approximation and is known to 
have limitations. For example, the Smoluchowski density 
decay amplitude in d = 1 differs from exact solutions by 
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a factor of 7r/2 13], and in multi-species reactions even 
the exponents can be wrong ^4| . 

However, Smoluchowski theory gives, surprisingly, the 
correct amplitude for the kA — > IA reaction at the upper 
critical dimension d c — 2/ (k — 1) for k = 2 and 3, as we 
demonstrate in this section. We first present the basic 
approach for two-particle reactions A + A — ► (0,j4) in 
some detail, as this is necessary for getting the various 
factors correct and for generalizing to the three-particle 
case. 



A. A + A -> (0, A) reaction in d = 2 

A coordinate system origin is attached to one of par- 
ticles. The motion of the other particles in this coordi- 
nate system remains diffusive with an effective diffusion 
constant D = 2D, as can be straightforwardly verified 
from the continuum Green's functions for the diffusion 
equation. The Smoluchowski theory is then built on the 
conditional density n(r,t), where r is the radial distance 
from the particle at the origin. First, n(r, f) is assumed 
to evolve via the diffusion equation, 



g-n(r,t) 



DV 2 n{r, t), 



(21) 



subject to boundary conditions n(R, t) = 0, where R is 
the radius of the fixed particle, and n(r — > oo, t) = n^, 
and initial condition n(r, 0) = n^. The solution of 
eq. (|21|) is used to determine the particle current j = 
— DVn(f,f), which in turn gives the flux of particles to- 
ward the origin, 



F = - $ <ISi<-j\,-,.- = --S.iD[j£ 



(22) 



r=R 



where Sd is the surface area of the d-dimensional unit 
hypersphere. The flux F then determines the decay rate 
of the probability p that the particle at the origin has not 
been visited by another particle, i.e., p = —Fp. 

This can be turned into an equation for the density 
decay by considering each particle in turn to be the fixed 
particle. Since each meeting of particles involves two 
particles and removes 2 — 1 particles, the density will 
decay according to 



dn 
~dt 



2-1. 



-Fn. 



(23) 



To close these equations, ««, is taken to be n(t). This 
amounts to a quasistatic approximation, since the time- 
dependence of the large r boundary condition is ne- 
glected when solving eq. J2U- The second approxima- 
tion in Smoluchowski theory is that the exterior parti- 
cles are treated as only diffusing, and reactions are re- 
incorporated in a mean-field way via the decaying density 
used as the large r boundary condition. 



Now consider the case of d = 2. For times t R 2 /D 
the solution to eq. (|21|l in the region R < r <C (Dt) 1 / 2 is 
given by 0] 



n(r, t) = 



lnt 



]n(r/R). 



(24) 



The resulting flux is F = 47rnooD/lni. Taking n c 
n(t) gives, via eq. (|2*3l . 



dn 
~dl 



4tt(2 - t)D 

; ' 

In* 



(25) 



where we have substituted back the lab frame diffusion 
constant. This has the form of the mean- field rate equa- 
tion with a time-dependent rate constant that decays as 
1/ In t. Eq. (|25|l results in an asymptotic density of the 
form n ~ A\nt/(Dt), with the amplitude matching the 
RG value, eq. 10), and the exact solution [(j. We believe 
this to be the first demonstration that Smoluchowski the- 
ory predicts the correct amplitude at the upper critical 
dimension. 

Finally, we comment on the case where the particles 
are not circular. The boundary condition for the con- 
ditional density remains n(r, t) = for points r on the 
particle surface. While we are unaware of a formal ex- 
act solution to eq. I|21|) for this case, we appeal to the 
quasistatic approach presented in [l(|- The conditional 
density is assumed to obey Laplace's equation in the re- 
gion exterior to the particle at the origin but within the 
diffusion radius Rd — a(Dt) 1 / 2 , where a is an arbitrary 
constant. The outer boundary condition is taken to be 
n(R]j, t) = rirx,. For a circular particle, this approach re- 
produces the known solution, eq. (|24|l . For a non-circular 
particle, it is straightforward to demonstrate via separa- 
tion of variables that the solution is unmodified apart 
from an additive constant (with respect to r) and ad- 
ditional subasymptotic terms that decay as powers of 
t. This suggests the asymptotic Smoluchowski flux and 
subsequent density are universal with respect to particle 
shape. 



B. 3^4 — > I A reaction in d — 1 

Two approaches have been used to extend Smolu- 
chowski theory to the case of three-particle, one species 
reactions in one dimension. Both approaches yielded den- 
sities with the logarithmic correction to the rate equation 
result. 

Krapivsky's method |10| is to consider a fixed parti- 
cle and then construct pseudo-particles in a d = 2 plane 
out of every possible pair of particles, i.e. real particles 
at points x\ and xi in the one-dimensional system would 
contribute a pseudo-particle at (x\ , x%). In this viewpoint 
a three-particle reaction corresponds to a pseudo-particle 
meeting the fixed particle at the origin. Krapivsky pro- 
poses to approximate the pseudo-particle dynamics as in- 
dependently diffusing particles, even though their motion 
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is correlated, with the net result that the problem reduces 
to the d = 2 two-particle reaction. The only difference is 
that the represents the pseudo-particle density, which 
is given by n(t) 2 , where n(t) is the d — 1 particle density. 
The corresponding equation becomes dtn ~ — Dn 3 / In t, 
which exhibits the expected n ~ y/\n.t/Dt behavior. 

The second approach, due to Oshanin, Stemmer, Lud- 
ing, and Blumen (OSLB) 3], is based on the three-point 
correlation function, which due to translational invari- 
ance is a function of two distances. OSLB show that 
this correlation function, in the absence of reactions, sat- 
isfies an anisotropic two-dimensional diffusion equation 
in the plane defined by the two arguments. Reactions 
are then incorporated into this pseudo-two-dimensional 
system via the Smoluchowski approach. 

We argue that at the level of the Smoluchowski ap- 
proximation these two approaches are the same. The 
three-point correlation function viewed in the 2D plane is, 
in fact, the average of Krapivsky's pseudoparticlcs over 
stochastic initial conditions and diffusion hops. Such an 
averaging is implicit in the Krapivsky method in going 
to the pseudoparticle diffusion equation, so at this point 
the two approaches should be formally identical. We now 
proceed to solve this system. 

First, we observe that Krapivsky's pseudo-particles in- 
deed satisfy the same anisotropic diffusion equation as 
the OSLB. Let j/i(t) = Xi(t) — Xi(0) represent the dis- 
placement at time t of the ith particle from its ini- 
tial position. The single particle Green's function is 
G(y, t) oc e~ y i( iDt \ where we neglect normalization fac- 
tors for clarity. We attach a coordinate origin to particle 
i = by introducing the variables m — y% — ya- Now 
consider a pair of exterior, i ^ particles, which we take 
to be i = 1 and 2. The combined Green's function in the 
Smoluchowski frame is given by 



G 2 (ui , u 2 , t) = J dy G(y Q , t) G{y x , t) G(y 2 , t) 



oc exp 



(26) 



The u\u 2 cross-term indicates anisotropic diffusion, 
which may be diagonalized by a 7r/4 rotation. Taking 
Vi = (iti + U2)/v2 and v 2 — {u\ — u 2 )/V2 gives 



G 2 (vi , V2, t) oc exp 



YIDt 



u 2 

ADt 



(27) 



By comparison with the single-particle Green's function 
we find that in the v\ direction the effective diffusion 
constant is D\ = 3D, while in the v 2 direction D 2 = D. 
This can be understood by noting that changes in yo, 
that is, the motion of the particle at the origin, affects 
Mi and i*2 identically, thus enhancing the motion of their 
sum relative to their difference. 

To make the Smoluchowski approximation one then 
solves this anisotropic diffusion equation in the v plane. 
We first rescale v\ = vi/y/3 and v 2 — v 2 to get isotropic 
diffusion. The rescaling affects the shape of the interior 
boundary condition, but as argued above, this should 



not affect the position-dependent part of the density at 
distances small compared to (Dt) 1 / 2 . Mapping back to 
the v plane we find 

n(vi,v 2 ,t) - (n 0O /lnt)ln(uJ/3 + t;2) + const. (28) 

in the region vi,v 2 <C (Dt) 1 / 2 . The particle current in 
the pseudo-particle plane is then 



(29) 
(30) 



(vivi + v 2 i) 2 ) 



2n x D 



Interestingly, current is directed radially inward despite 
the anisotropy. The flux through a circular region en- 
compassing the fixed particle is 



F = 



2nooD_ 
hit 



dB 



i cos 2 6 + sin 2 9 



4v / 3~7m 0O Z) 
lot 



(31) 



Now we generate an equation for the density n(t) of 
our one-dimensional system. As before, we generate an 
equation for the density by considering each particle in 
turn to be fixed at the origin. This brings a factor of 
(3 — £)/3 since each reaction removes 3 — I particles and 
is counted three times. Furthermore, since each particle 
pair contributes two pseudo-particles, i.e., at (xi, x 2 ) and 
(x2,xi), the flux of pseudoparticles to the origin double 
counts the reactions with the particle at the origin. The 
net result is 



dn 
~dt 



-Fn 



2^3(3 - £)ttD 



31ni 



(32) 



where in the last step we have taken rioo — n(t) 2 . As 
with the A + A — > (0, -A) reactions in d = 2, this has the 
rate equation form with a l/ln£ reaction rate. The re- 
sulting density is of the expected form n ~ Ay/]nt/ (Dt), 
with the amplitude A matching the RG result, eq. J2J- 
However, we note that the universal leading corrections 
predicted by the RG are absent, i.e., according Smolu- 
chowski theory — 0. 

Our result for the Smoluchowski amplitudes differ 
quantitatively from Krapivsky |10| , who did not consider 
the anisotropic nature of the diffusion, and from OSLB 
[3| , who omitted the factor of three due to the number of 
particles removed in 3 A — > and the factor of 1/2 due to 
double-counting the rate of reactions at the origin. With 
these factors included, the OSLB result becomes equiva- 
lent to eq. (|32"|) . 



VI. SUMMARY 

The present work was motivated by the apparent dis- 
crepancy between the simulations and RG predictions. 
This discrepancy was resolved by demonstrating that the 
densities measured by simulation, which appear to be 
asymptotic in plots such as Fig. ^ are in fact consistent 
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with the RG predictions of significant subasymptotic con- 
tributions, as demonstrated in Figs. [3 El and^] Unfor- 
tunately, we also conclude that a more precise test of the 
RG predictions by direct comparison to simulations is 
not possible, owing to the slow decay of the higher order 
terms in a 1/ \/hit expansion. 

In view of this, it is significant that exact solutions are 
available for the A + A — > (0, A) reactions in d = d c = 2, 
which do provide a strong test of the RG results. It is also 
significant that relations predicted via the field rescaling 
transformation between pure and mixed reactions can be 
tested by simulations and appear to hold with no slow 
transients. Given that these k = 2 amplitudes are cor- 
rect, that the field rescaling transformation predictions 
hold, that there is no discrepancy with simulation in the 
pure 3^4 — > IA reactions, and that in principle no approx- 
imation is made in the RG treatment, it seems reasonable 
to believe the RG amplitudes reported here are exact re- 
sults. 

As is the usual benefit of an RG calculation, these am- 
plitudes are also demonstrated to be universal, thus the 
predictions are subject to testing by future exact solu- 



tions for any particular realization of the dynamics. 

A final result of our study is the demonstration that 
there is a single Smoluchowski theory for the 3^4 — > IA 
reaction in one dimension, and that this approximate the- 
ory indeed yields the exact asymptotic density. Since 
this property is shared by the A + A — > (0, A) reaction in 
d = 2, it appears to be a general feature of Smoluchowski 
theory that it succeeds quantitatively at the upper criti- 
cal dimension. It may be possible to find an underlying 
explanation for this property, which would be an inter- 
esting direction for future work. 
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